Antenna Coupling Component Measurement Tool Having a Rotating Antenna Configuration

ABSTRACT

Disclosed herein are electromagnetic resistivity logging systems and methods that employ an antenna configuration having at most two transmitter or receiver antenna orientations that rotate relative to the borehole. The measurements made by this reduced-complexity antenna configuration enable the determination of at least seven components of a coupling matrix, which may be determined using a linear system of equations that express the azimuthal dependence of the measurements. For increased reliability, measurement averaging may be performed in azimuthally spaced bins. The coupling matrix components can then be used as the basis for determining logs of various formation parameters, including vertical resistivity and anisotropy.

CROSS REFERENCE TO RELATED APPLICATIONS

The present application is a divisional application of parent U.S. patent application Ser. No. 12/294,557, filed Sep. 25, 2008, by Michael S. Bittar and Guoyu (David) Hu, entitled “Antenna Coupling Component Measurement Tool Having a Rotating Antenna Configuration, which claims priority to PCT Application No. PCT/US2006/062149, filed Dec. 15, 2006, entitled “Antenna Coupling Component Measurement Tool Having a Rotating Antenna Configuration”, by the same Inventors.

BACKGROUND

The basic principles and techniques for electromagnetic logging for earth formations are well known. For example, induction logging to determine the resistivity (or its inverse, conductivity) of earth formations adjacent a borehole has long been a standard and important technique in the search for and recovery of subterranean petroleum deposits. In brief, a transmitter transmits an electromagnetic signal that passes through formation materials around the borehole and induces a signal in ore or more receivers. The amplitude and/or phase of the receiver signals are influenced by the formation resistivity, enabling resistivity measurements to be made. The measured signal characteristics and/or formation properties calculated therefrom are recorded as a function of the tool's depth or position in the borehole, yielding a formation log that can be used by analysts.

Note, however, that the resistivity of a given formation may be isotropic (equal in all directions) or anisotropic (unequal in different directions). In electrically anisotropic formations, the anisotropy is generally attributable to extremely fine layering during the sedimentary build-up of the formation. Hence, in a formation coordinate system oriented such that the x-y plane is parallel to the formation layers and the z axis is perpendicular to the formation layers, resistivities R_(X) and R_(Y) in directions x and y, respectively, tend to be the same, but resistivity R_(Z) in the z direction is different. Thus, the resistivity in a direction parallel to the plane of the formation (i.e., the x-y plane) is known as the horizontal resistivity, R_(H), and the resistivity in the direction perpendicular to the plane of the formation (i.e., the z direction) is known as the vertical resistivity, Rv. The index of anisotropy, η, is defined as η=[R_(V)/R_(H)]^(1/2).

As a further complication to measuring formation resistivity, boreholes are generally not perpendicular to formation beds. The angle between the axis of the well bore and the orientation of the formation beds (as represented by the normal vector) has two components. These components are the dip angle and the strike angle. The dip angle is the angle between the borehole axis and the normal vector for the formation bed. The strike angle is the direction in which the boreholes axis “leans away from” the normal vector. (These will be defined more rigorously in the detailed description.)

Electromagnetic resistivity logging measurements are a complex function of formation resistivity, formation anisotropy, and the formation dip and strike angles, which may all be unknown. Logging tools that fail to account for one or more of these parameters may provide measurement quality that is less than ideal. Conversely, tools that can be used to measure each of these parameters will provide improved resistivity measurements. Moreover, tools that are able to provide dip and strike measurements along with azimuthal orientation information, can be used for geosteering. (Geosteering is a process in which drill engineers adjust the drilling direction to increase the borehole's exposure to a hydrocarbon-bearing formation (the “payzone”).)

BRIEF DESCRIPTION OF THE DRAWINGS

A better understanding of the various disclosed embodiments can be obtained when the following detailed description is considered in conjunction with the following drawings, in which:

FIG. 1 shows an illustrative logging while drilling environment including dipping formation beds;

FIG. 2 shows an illustrative wireline logging environment including dipping formation beds;

FIG. 3 shows a relationship between the orientation of a borehole and a dipping formation bed;

FIG. 4 shows a hypothetical antenna arrangement for a tool having an orthogonal triaxial transmitter and two orthogonal triaxial receivers;

FIG. 5 shows angles for defining the orientation of a tilted antenna;

FIG. 6 is a block diagram of an illustrative electronics module for an electromagnetic resistivity tool;

FIG. 7 shows an illustrative electromagnetic resistivity logging tool having tilted transmitter and receiver antennas;

FIG. 8 is a flow diagram of an illustrative electromagnetic resistivity logging method;

FIG. 9 shows an illustrative electromagnetic resistivity logging tool having parallel tilted transmitter and receiver antennas;

FIG. 10 shows an illustrative electromagnetic resistivity logging tool having transmitters tilted at a first orientation and receiver antennas tilted at a second orientation;

FIG. 11 shows an illustrative electromagnetic resistivity logging tool having both parallel and non-parallel tilted transmitter and receiver antennas;

FIG. 12 shows an illustrative electromagnetic resistivity logging tool having co-located tilted receiver antennas;

FIG. 13 shows the division of a borehore circumference into azimuthal bins; and

FIG. 14 shows an illustrative electromagnetic resistivity logging tool having compensated measurements.

While the invention is susceptible to various modifications and alternative forms, specific embodiments thereof are shown by way of example in the drawings and will herein be described in detail. It should be understood, however, that the drawings and detailed description thereto are not intended to limit the invention to the particular form disclosed, but on the contrary, the intention is to cover all modifications, equivalents and alternatives falling within the spirit and scope of the appended claims.

DETAILED DESCRIPTION

Disclosed herein are electromagnetic resistivity logging systems and methods that employ rotation of an azimuthally sensitive antenna configuration to obtain more elements of a coupling matrix than could be obtained from the static antenna configuration alone. Such exploitation of the rotation enables a complete antenna coupling matrix to be obtained with a tool having only two antenna orientations. With the complete antenna coupling matrix, the various resistivity measurement parameters can be measured, including formation resistivity, formation anisotropy, and formation dip and strike angles. Moreover, when the coupling matrix values are combined with orientation information from the tool, geosteering or boundary detection signals may be derived. Such benefits may be obtained with a reduced cost and greater reliability due to the reduced number of antennas.

The disclosed tool configurations and operations are best understood in the context of the larger systems in which they operate. Accordingly, an illustrative logging while drilling (LWD) environment is shown in FIG. 1. A drilling platform 2 supports a derrick 4 having a traveling block 6 for raising and lowering a drill string 8. A kelly 10 supports the drill string 8 as it is lowered through a rotary table 12. A drill bit 14 is driven by a downhole motor and/or rotation of the drill string 8. As bit 14 rotates, it creates a borehole 16 that passes through various formations 18. A pump 20 circulates drilling fluid through a feed pipe 22 to kelly 10, downhole through the interior of drill string 8, through orifices in drill bit 14, back to the surface via the annulus around drill string 8, and into a retention pit 24. The drilling fluid transports cuttings from the borehole into the pit 24 and aids in maintaining the borehole integrity.

An electromagnetic resistivity logging tool 26 is integrated into the bottom-hole assembly near the bit 14. As the bit extends the borehole through the formations, logging tool 26 collects measurements relating to various formation properties as well as the tool orientation and position and various other drilling conditions. (The orientation measurements may be performed using an azimuthal orientation indicator, which may include magnetometers, inclinometers, and/or accelerometers, though other sensor types such as gyroscopes may be used. In some embodiments, the tool includes a 3-axis fluxgate magnetometer and a 3-axis accelerometer.) The logging tool 26 may take the form of a drill collar, i.e., a thick-walled tubular that provides weight and rigidity to aid the drilling process. A telemetry sub 28 may be included to transfer tool measurements to a surface receiver 30 and to receive commands from the surface receiver.

More preferably, rotational position indicator 70 may contain both a 3-axis fluxgate magnetometer and a 3-axis accelerometer. As is known in the art, the combination of those two sensor systems enables the measurement of the toolface, inclination, and azimuth orientation angles of the borehole. The toolface and hole inclination angles are calculated from the accelerometer sensor output. The magnetometer sensor outputs are used to calculate the hole azimuth. With the toolface, the hole inclination, and the hole azimuth information, a tool in accordance with the present invention can be used to steer the bit to the desirable bed. Specifically, the response difference or the response ratio can be used effectively to enter a desired payzone or to stay within the payzone of interest.

At various times during the drilling process, the drill string 8 may be removed from the borehole as shown in FIG. 2. Once the drill string has been removed, logging operations can be conducted using a wireline logging tool 34, i.e., a sensing instrument sonde suspended by a cable 42 having conductors for transporting power to the tool and telemetry from the tool to the surface. A resistivity imaging portion of the logging tool 34 may have centralizing arms 36 that center the tool within the borehole as the tool is pulled uphole. A logging facility 44 collects measurements from the logging tool 34, and includes computing facilities for processing and storing the measurements gathered by the logging tool.

FIG. 1 shows that the formations 18 are not perpendicular to the borehole, which may occur naturally or due to directional drilling operations. The borehole may have a coordinate system 50 defined in accordance with the borehole's long axis (the z axis) and the north side (or alternatively, the high side) of the hole (the x-axis). The formations 18, when characterized as a plane, may have a coordinate system 51 defined in accordance with the normal to the plane (the z″ axis) and the direction of steepest descent (the x″-axis). As shown in FIG. 3, the two coordinate systems are related by two rotations. Beginning with the borehole's coordinate system (x,y,z), a first rotation of angle γ is made about the z axis. The resulting coordinate system is denoted (x′,y′,z′). Angle γ is the relative strike angle, which indicates the direction of the formation dip relative to the borehole's coordinate system. A second rotation of angle α is then made about the y′ axis. This aligns the borehole coordinate system with the formation coordinate system. Angle α is the relative dip angle, which is the slope angle of the beds relative to the long axis of the borehole.

The vertical resistivity is generally found to be the resistivity as measured perpendicular to the plane of the formation, and the horizontal resistivity is the resistivity as measured within the plane of the formation. Determination of each of these parameters (dip angle, strike angle, vertical resistivity, and horizontal resistivity) is desirable.

FIG. 4 shows a hypothetical antenna configuration for a multi-component electromagnetic resistivity logging tool. (The electromagnetic resistivity logging tool may be embodied as a wireline tool and as a logging while drilling tool.) A triad of transmitter coils T_(X), T_(Y) and T_(Z), each oriented along a respective axis, is provided. At least one triad of similarly oriented receiver coils R_(1X), R_(1Y), and R_(1Z) is also provided. For receive signal measurements relative to the amplitude and phase of the transmit signal (sometimes called “absolute” measurements), only one receiver triad would be used. A second triad of similarly oriented receiver coils pairs R_(2X), R_(2Y), and R_(2Z) may also provided when differential measurements are desired (e.g., a signal amplitude ratio or a phase difference between receiver coils oriented along a given axis). Differential measurements may offer increased spatial resolution.

Moran and Gianzero, in “Effects of Formation Anisotropy on Resistivity Logging Measurements” Geophysics, Vol. 44, No. 7, p. 1266 (1979), noted that the magnetic field h in the receiver coils can be represented in terms of the magnetic moments m at the transmitters and a coupling matrix C:

h=Cm  (1)

In express form, equation (1) is:

$\begin{matrix} {{\begin{bmatrix} H_{x} \\ H_{y} \\ H_{z} \end{bmatrix} = {\begin{bmatrix} C_{xx} & C_{xy} & C_{xz} \\ C_{yx} & C_{yy} & C_{zz} \\ C_{zx} & C_{zy} & C_{zz} \end{bmatrix}\begin{bmatrix} M_{x} \\ M_{y} \\ M_{z} \end{bmatrix}}},} & (2) \end{matrix}$

where M_(X), M_(Y), and M_(Z) are the magnetic moments (proportional to transmit signal strength) created by transmitters T_(X), T_(Y), and T_(Z), respectively. H_(X), H_(Y), H_(Z) are the magnetic fields (proportional to receive signal strength) at the receiver antennas R_(X), R_(Y), and R_(Z), respectively.

In the antenna configuration of FIG. 4, if each transmitter is fired in turn, and signal measurements are made at each receiver in response to each firing, nine absolute or differential measurements are obtained. These nine measurements enable the determination of a complete coupling matrix C. (C_(IJ)=a_(IJ)V_(I) ^(J), where I is the index for receiver R_(X), R_(Y), or R_(Z), J is the index for transmitter T_(X), T_(Y), or T_(Z), a_(IJ) is a constant determined by the tool design, and V_(I) ^(J) is a complex value representing the signal amplitude and phase shift measured by receiver I in response to the firing of transmitter J.) Knowledge of the complete coupling matrix enables the determination of dip angle, strike angle, vertical resistivity, and horizontal resistivity. A number of techniques may be used to determine these parameters. For example, dip and strike angle may be determined from coupling matrix values as explained by Li Gao and Stanley Gianzero, U.S. Pat. No. 6,727,706 “Virtual Steering of Induction Tool for Determination of Formation Dip Angle”. Given these angles, vertical and horizontal resistivity can be determined in accordance with equations provided by Michael Bittar, U.S. Pat. No. 7,019,528 “Electromagnetic Wave Resistivity Tool Having a Tilted Antenna for Geosteering within a Desired Payzone”. Alternatively, a simultaneous solution for these parameters may be found as described in the Bittar reference.

FIG. 5 shows two angles that may be used to specify the orientation of a coil antenna. The coil antenna may be considered as residing in a plane having a normal vector. Tilt angle θ is the angle between the longitudinal axis of the tool and the normal vector. Azimuth angle β is the angle between the projection of the normal vector in the X-Y plane and the tool scribe line. Alternatively, in the downhole context, azimuthal angle β may represent the angle between projection of the normal vector in the X-Y plane and the x-axis of the borehole coordinate system.

It is noted that three transmitter antenna orientations and three receiver antenna orientations are employed in the antenna configuration of FIG. 4. It has been discovered that when tool rotation is exploited, it is possible to determine the full coupling matrix with only one transmitter and two receiver antenna orientations (or equivalently, one receiver and two transmitter antenna orientations). Moreover, with certain assumptions about the configuration of the formation, one transmitter and receiver antenna orientation may be sufficient.

Before considering various tools having specific antenna configurations, the electronics common to each tool are described. FIG. 6 shows a functional block diagram of the electronics for a resistivity tool. The electronics include a control module 602 that is coupled to an analog switch 604. Analog switch 604 is configured to drive any one of the transmitter coils T₁, T₂, T₃, T₄ with an alternating current (AC) signal from a signal source 606. In some embodiments, the signal source provides radio frequency signals. The control module 602 preferably selects a transmitter coil, pauses long enough for transients to die out, then signals data storage/transmit module 610 to accept an amplitude and phase sample of the signals received by each of the receivers. The control module 602 preferably repeats this process sequentially for each of the transmitters. The amplitude and phase shift values are provided by amplitude and phase shift detector 608 which is coupled to each of the receiver coils R₁ and R₂ for this purpose.

Control module 602 may process the amplitude and phase shift measurements to obtain compensated measurements and/or measurement averages. The raw, compensated, or averaged measurements, may be transmitted to the surface for processing to determine dip and strike angles, vertical and horizontal resistivity, and other information such as (i) distance to nearest bed boundary, (ii) direction of nearest bed boundary, and (iii) resistivity of any nearby adjacent beds. The data storage/transmitter module 610 may be coupled to telemetry unit 28 (FIG. 1) to transmit signal measurements to the surface. Telemetry unit 28 can use any of several known techniques for transmitting information to the surface, including but not limited to (1) mud pressure pulse; (2) hard-wire connection; (3) acoustic wave; and (4) electromagnetic waves.

FIG. 7 shows an electromagnetic resistivity logging tool 702 having only two receiver antenna orientations. The tool 702 is provided with one or more regions 706 of reduced diameter. A wire coil 704 is placed in the region 706 and in some embodiments is spaced away from the surface of subassembly 702 by a constant distance. To mechanically support and protect the coil 704, a non-conductive filler material (not shown) such as epoxy, rubber, or ceramic may be used in the reduced diameter regions 706. Coil 704 is a transmitter coil, and coils 710 and 712 are receiving coils. In operation, transmitter coil 704 transmits an interrogating electromagnetic signal which propagates through the borehole and surrounding formation. Receiver coils 710, 712 detect the interrogating electromagnetic signal and provide a measure of the electromagnetic signal's amplitude attenuation and phase shift. For differential measurements additional receiver coils parallel to coils 710, 712 may be provided at an axially-spaced distance. From the absolute or differential amplitude attenuation and phase shift measurements, the coupling matrix components can be determined and used as the basis for determining formation parameters and as the basis for geosteering.

The transmitter coil 704 may be spaced approximately 30 inches from the receiver coils 710, 712. The transmitter and receiver coils may comprise as little as one loop of wire, although more loops may provide additional signal power. The distance between the coils and the tool surface is preferably in the range from 1/16 inch to ¾ inch, but may be larger. Transmitter coil 704 and receiver coil 712 may each have a tilt angle of about 45° and aligned with the same azimuth angle, while receiver coil 710 may have a tilt angle of about 45° and an azimuth 180° apart from receiver coil 712 (or equivalently, a tilt angle of minus 45° at the same azimuth angle as receiver coil 712).

The signal measured by a tilted receiver in response to the firing of a tilted transmitter can be expressed in terms of the signals V_(I) ^(J) that would be measured by the tool of FIG. 4. When both transmitter and receiver coils are oriented at the same azimuth angle β, the tilted receiver signal V_(R) is

$\begin{matrix} {V_{R} = {{\begin{bmatrix} {\sin \mspace{14mu} \theta_{T}\mspace{14mu} \cos \mspace{14mu} \beta} \\ {\sin \mspace{14mu} \theta_{T}\mspace{14mu} \sin \mspace{14mu} \beta} \\ {\cos \mspace{14mu} \theta_{T}} \end{bmatrix}\begin{bmatrix} V_{x}^{x} & V_{x}^{y} & V_{x}^{z} \\ V_{y}^{x} & V_{y}^{y} & V_{y}^{z} \\ V_{z}^{x} & V_{z}^{y} & V_{z}^{z} \end{bmatrix}}\begin{bmatrix} {\sin \mspace{14mu} \theta_{R}\mspace{14mu} \cos \mspace{14mu} \beta} \\ {\sin \mspace{14mu} \theta_{R}\mspace{14mu} \sin \mspace{14mu} \beta} \\ {\cos \mspace{14mu} \theta_{R}} \end{bmatrix}}} & (3) \end{matrix}$

where θ_(T) is the tilt angle of the transmitter and θ_(R) is the tilt angle of the receiver. In written-out form, the tilted receiver angle is:

$\begin{matrix} {V_{R} = \left\lbrack {{V_{x}^{x}\mspace{14mu} \sin \mspace{14mu} \theta_{T}\mspace{14mu} \sin \mspace{14mu} \theta_{R}\mspace{14mu} \cos^{2}\mspace{14mu} \beta} + {V_{y}^{x}\mspace{14mu} \sin \mspace{14mu} \theta_{T}\mspace{14mu} \sin \mspace{14mu} \theta_{R}\mspace{14mu} \sin \mspace{14mu} \beta \mspace{14mu} \cos \mspace{14mu} \beta} + {V_{z}^{x}\mspace{14mu} \cos \mspace{14mu} \theta_{T}\mspace{14mu} \sin \mspace{14mu} \theta_{R}\mspace{14mu} \cos \mspace{14mu} \beta} + {V_{x}^{y}\mspace{14mu} \sin \mspace{14mu} \theta_{T}\mspace{14mu} \sin \mspace{14mu} \theta_{R}\mspace{14mu} \sin \mspace{14mu} \beta \mspace{14mu} \cos \mspace{14mu} \beta} + {V_{y}^{y}\mspace{14mu} \sin \mspace{14mu} \theta_{T}\mspace{14mu} \sin \mspace{14mu} \theta_{R}\mspace{14mu} \sin^{2}\mspace{14mu} \beta} + {V_{x}^{y}\mspace{14mu} \cos \mspace{14mu} \theta_{T}\mspace{14mu} \sin \mspace{14mu} \theta_{R}\mspace{14mu} \sin \mspace{14mu} \beta} + {V_{x}^{z}\mspace{14mu} \sin \mspace{14mu} \theta_{T}\mspace{14mu} \cos \mspace{14mu} \theta_{R}\mspace{14mu} \cos \mspace{14mu} \beta} + {V_{y}^{z}\mspace{14mu} \sin \mspace{14mu} \theta_{T}\mspace{14mu} \cos \mspace{14mu} \theta_{R}\mspace{14mu} \sin \mspace{14mu} \beta} + {V_{z}^{z}\mspace{14mu} \cos \mspace{14mu} \theta_{T}\mspace{14mu} \cos \mspace{14mu} \theta_{R}}} \right\rbrack} & (4) \end{matrix}$

It is noted that the azimuth angle dependence for V_(y) ^(x) and V_(x) ^(y) is the same (sin β cos β), the dependence for V_(z) ^(x) and V_(x) ^(z) is the same (cos β), and the dependence for V_(z) ^(y) and V_(y) ^(z) is the same (sin β). It may be assumed that cross coupling components V_(y) ^(x) and V_(x) ^(y) are equal, but such an assumption is not desirable for the remaining cross components (at least not in dipping anisotropic beds). In that situation, the cross coupling components cannot be independently determined from a single rotating tilted transmitter-receiver pair such as transmitter coil 704 and receiver coil 712. (Note however that the diagonal elements can still be calculated.) A second transmitter or receiver coil (e.g., receiver coil 712) may be employed to provide an independent set of equations that enable the cross-coupling values to be determined.

A linear system of equations for the measurements made by rotating tool 702 is provided below. In deriving these equations, the borehole coordinate system is chosen so that the x-axis aligns with the dip azimuth of the surrounding formation, causing V_(y) ^(x) and V_(x) ^(y) to be zero. (The dip azimuth may be found by determining the azimuth angles at which the imaginary component of the receive signal reaches its minimum magnitude. Alternatively, these components can be assumed to be equal and can be retained in the linear system of equations.) The signals measured by receivers R1 (coil 712) and R2 (coil 710) at various azimuthal angles β₁-β_(N) are:

$\begin{matrix} {\; {{{V_{R\; 1}\left( \beta_{1} \right)} = {\frac{1}{2}\left\lbrack {{V_{x}^{x}\mspace{14mu} \cos^{2}\mspace{14mu} \beta_{1}} + {V_{z}^{x}\mspace{14mu} \cos \mspace{14mu} \beta_{1}} + {V_{y}^{y}\mspace{14mu} \sin^{2}\mspace{14mu} \beta_{1}} + {V_{z}^{y}\mspace{14mu} \sin \mspace{14mu} \beta_{1}} + {V_{y}^{z}\mspace{14mu} \sin \mspace{14mu} \beta_{1}} + V_{z}^{z}} \right\rbrack}},}} & (5.1) \\ {{{V_{R\; 1}\left( \beta_{2} \right)} = {\frac{1}{2}\left\lbrack {{V_{x}^{x}\mspace{14mu} \cos^{2}\mspace{14mu} \beta_{2}} + {V_{z}^{x}\mspace{14mu} \cos \mspace{14mu} \beta_{2}} + {V_{y}^{y}\mspace{14mu} \sin^{2}\mspace{14mu} \beta_{2}} + {V_{z}^{y}\mspace{14mu} \sin \mspace{14mu} \beta_{2}} + {V_{x}^{z}\mspace{14mu} \cos \mspace{14mu} \beta_{2}} + {V_{y}^{z}\mspace{14mu} \sin \mspace{14mu} \beta_{2}} + V_{z}^{z}} \right\rbrack}},\ldots} & (5.2) \\ {{{V_{R\; 1}\left( \beta_{N} \right)} = {\frac{1}{2}\left\lbrack {{V_{x}^{x}\mspace{14mu} \cos^{2}\mspace{14mu} \beta_{N}} + {V_{z}^{x}\mspace{14mu} \cos \mspace{14mu} \beta_{N}} + {V_{y}^{y}\mspace{14mu} \sin^{2}\mspace{14mu} \beta_{N}} + {V_{z}^{y}\mspace{14mu} \sin \mspace{14mu} \beta_{N}} + {V_{x}^{z}\mspace{14mu} \cos \mspace{14mu} \beta_{N}} + {V_{y}^{z}\mspace{14mu} \sin \mspace{14mu} \beta_{N}} + V_{z}^{z}} \right\rbrack}},} & \left( {5.N} \right) \\ {{{V_{R\; 2}\left( \beta_{1} \right)} = {\frac{1}{2}\left\lbrack {{{- V_{x}^{x}}\mspace{14mu} \cos^{2}\mspace{14mu} \beta_{1}} - {V_{z}^{x}\mspace{14mu} \cos \mspace{14mu} \beta_{1}} - {V_{y}^{y}\mspace{14mu} \sin^{2}\mspace{14mu} \beta_{1}} - {V_{z}^{y}\mspace{14mu} \sin \mspace{14mu} \beta_{1}} + {V_{x}^{z}\mspace{14mu} \cos \mspace{14mu} \beta_{1}} + {V_{y}^{z}\mspace{14mu} \sin \mspace{14mu} \beta_{1}} + V_{z}^{z}} \right\rbrack}},\ldots} & \left( {{5.N} + 1} \right) \\ {{V_{R\; 2}\left( \beta_{1} \right)} = {{\frac{1}{2}\left\lbrack {{{- V_{x}^{x}}\mspace{14mu} \cos^{2}\mspace{14mu} \beta_{N}} - {V_{z}^{x}\mspace{14mu} \cos \mspace{14mu} \beta_{N}} - {V_{y}^{y}\mspace{14mu} \sin^{2}\mspace{14mu} \beta_{N}} - {V_{z}^{y}\mspace{14mu} \sin \mspace{14mu} \beta_{N}} + {V_{x}^{z}\mspace{14mu} \cos \mspace{14mu} \beta_{N}} + {V_{y}^{z}\mspace{14mu} \sin \mspace{14mu} \beta_{N}} + V_{z}^{z}} \right\rbrack}.}} & \left( {5.2N} \right) \end{matrix}$

Thus the linear system provided above has seven unknowns and 2N equations. It is expected that the tool would perform measurements for at least 10 different azimuthal angles, and may perform measurements for 16 or 32 evenly-divided azimuthal bins as shown in FIG. 13. The orthogonal voltage components V_(I) ^(J) can be readily determined using a linear least squares algorithm.

FIG. 8 shows an illustrative logging method which may be performed by a downhole controller, by a surface computing facility that receives measurements from the tool, or performed cooperatively by both. In block 802 an initial transmitter is selected. (Multi-transmitter tools are discussed further below.) In block 804, the selected transmitter is fired, and the amplitude and phase of each receiver's response is measured. The tool's position and orientation are also captured and used to associate the receiver response measurements with an azimuthal bin. (An azimuthal bin has both an angular extent and an axial extent.) In block 806, the current measurements are used to update the average response for each receiver for the given bin.

In block 808, a test is made to determine whether additional measurements are needed or will be forthcoming at the current borehole position. For example, in tools having multiple transmitters, it is desired to have measurements from each transmitter. Other reasons for needing additional measurements include having a desired number of measurements within each azimuthal bin before additional processing is performed, or having at least a given number of azimuthally different measurements before additional processing is performed. If additional measurements at the current position are expected, the additional processing may be postponed until all the relevant measurements have been collected.

Once a sufficient number of measurements have been obtained at a given position in the borehole, the method continues with optional block 810. In block 810, compensated measurements are calculated. Compensation is most suitable for differential receiver configurations like that shown in FIG. 14. Tool 1402 includes a pair of receivers 1410, 1412 located between equally spaced transmitters 1408, 1414. In response to the firing of the first transmitter 1408, receivers 1410 and 1412 detect attenuation and phase shift values A₁,Φ₁ and A₂,Φ₂, from which differential measurements can be determined (e.g., (Φ₂-Φ₁), or (logA₁-logA₂)). In response to the firing of the second transmitter 1414, receivers 1410 and 1412 detect attenuation and phase shift values A₄,Φ₄ and A₃,Φ₃, from which differential measurements can be determined (e.g., (Φ₄-Φ₃)). In each azimuthal bin, the differential measurement responses to the opposing transmitters can then be averaged together to obtain a compensated measurement, i.e., a measurement in which fixed biases in the electronics are canceled.

In block 812, the orthogonal antenna couplings are calculated based on the compensated measurements or on the average differential or average absolute measurements of attenuation and phase shift in each azimuthal bin. A least squares solution to a linear system of equations (such as that provided in equations 5.1 through 5.2N above) is calculated to find these coupling values. In block 814, the formation parameters of interest are calculated based on the orthogonal antenna couplings. The calculation of formation parameters may employ simultaneous inversion, or alternatively, some of the formation parameters may be fixed using outside information, thereby simplifying the solution process for the remaining parameters.

In optional block 816, a real-time log of one or more formation parameters (e.g., horizontal resistivity and anisotropy) is updated with the newly calculated parameter values from block 814. The log associates the calculated values with a depth or axial position within the borehole. Optionally, the information may also be associated with azimuthal orientations to generate a borehole image of azimuthal resistivity.

In block 818 a check is made to determine if logging information is available (or will become available) for additional positions within the borehole. If so, the process begins again with block 802. Otherwise, the process terminates.

FIGS. 7, 9-12 and 14 show various antenna configurations for illustrative electromagnetic resistivity logging tools that may be suitable alternatives for a multi-component logging tool such as that described with respect to FIG. 4. In FIG. 7, the tool includes a single tilted transmitter antenna 704 and at least one tilted receiver antenna 712 that is parallel to the transmitter antenna. The use of only a single tilted transmitter and receiver antenna orientation enables the determination of the diagonal coupling matrix components, and under assumptions of cross-component equality, enables the determination of the cross-coupling components too. However, in a preferred embodiment, a second tilted receiver antenna 710 is provided to enable determination of cross-coupling components while only assuming that the XY and YX components are equal. Thus tools having a single transmitter and two receiver antenna orientations (or a single receiver and two transmitter antenna orientations) may provide the best tradeoff between complexity and completeness.

Tool 702 may have several variations. In a first variation, the roles of transmitter and receiver are exchanged, so that transmitters coils 710 and 712 are alternately fired and the response of receiver coil 704 is measured. As additional variations, the crossed antenna coils 710 and 712 may be equally spaced in opposite directions from the antenna coil 704. Antenna coils 710 and 712 may retain their role as receivers responding to signals from transmitter coil 704, or the roles may again be exchanged.

Tool 702 is intended for absolute measurements (i.e., attenuation and phase shift are measured relative to the transmit signal). FIG. 9 shows an illustrative tool 902 intended for compensated differential measurements with different transmitter-receiver spacings. Tilted receiver coils 910 and 912 are parallel, and may have a tilt angle of about 45° and a spacing of about 8 inches. A first pair of parallel transmitters coils 908 and 914 are equally spaced from the receiver coil midpoint by about 32 inches, and are shown oriented parallel to the receiver antennas. In an alternative embodiment, the parallel transmitter coils 908 and 914 may have any tilt angle, including a tilt angle of zero (co-axial orientation). A second pair of transmitter coils 904 and 916 are equally spaced from the receiver coil midpoint by about 48 inches, and are parallel to the first pair of transmitter coils. The greater transmitter-receiver spacing enables the electromagnetic signals to provide measurements at greater penetration depths, enabling more accurate formation resistivity measurements. Since all of the transmitter coils are parallel and the receiver coils are parallel, tool 902 offers only a single transmitter antenna orientation and a single receiver orientation, meaning that the full coupling matrix can be determined only with a presumption of cross-coupling equalities. In one variation, tool 902 employs only a single receiver to absolute rather than differential measurements.

FIG. 10 shows an illustrative tool 1002 which also has parallel transmitter coils and parallel receiver coils in an antenna arrangement intended for making compensated differential measurements. The parallel receiver coils 1010 and 1012 are shown having a tilt of about 45° in an azimuth opposite to that of the transmitter coils. However, since only a single transmitter antenna orientation and a single receiver antenna orientation are provided, tool 1002 can enable determination of the full coupling matrix only with a presumption of cross-coupling equalities. In one uncompensated variation, tool 1002 omits transmitter coils 914 and 916.

FIG. 11 shows an illustrative tool 1102 having parallel transmitter coils, but having receiver coils 1110 and 1112 tilted about 45° in opposite azimuthal directions. Since tool 1102 provides two receiver orientations, determination of the full coupling matrix is possible. The configuration of tool 1102 enables compensated differential measurements to be performed. In an uncompensated variation, transmitters coils 914 and 916 are omitted.

FIG. 12 shows an illustrative tool 1202 having two transmitter antenna orientations and two receiver antenna orientations. Receiver antennas 1210 and 1212 are co-located and tilted about 45° in opposite azimuthal directions. A first transmitter coil pair includes transmitter coils 1208 and 914 equally spaced in opposite axial directions from the receiver antennas by about 32 inches, and a second transmitter coil pair includes coils 1204 and 916 equally spaced in opposite axial directions from the receiver antennas by about 48 inches. Each transmitter coil pair has transmitter coils tilted in opposite azimuthal directions. The configuration of FIG. 12 may be employed to make compensated differential measurements. In a non-differential variation, one of the receiver coils is omitted. In a second variation, the transmitter and receiver roles are exchanged, possibly with the omission of coil 1212. The role exchange between transmitter and receiver coils may be performed for each of the illustrative tools disclosed herein.

FIG. 14 shows an illustrative tool 1402 having a pair of parallel tilted receiver antennas 1410 and 1412, with a pair of parallel co-axial transmitter antennas 1402 and 1414 equally spaced from the midpoint of the receiver antennas. Though this and other disclosed antenna configurations have only a single transmitter orientation and a single receiver orientation, they are nevertheless sensitive to anisotropy in dipping beds and can be used to calculate the anisotropy.

As previously mentioned, components of the coupling matrix C may be used as a basis for geosteering. With a properly oriented borehole coordinate system, the z-axis is indicative of the direction to the bed boundary, and the C_(XZ) and C_(ZX) components are useful for determining the boundary's proximity.

Numerous variations and modifications will become apparent to those skilled in the art once the above disclosure is fully appreciated. For example, the foregoing disclosure describes numerous antenna configurations in the context of a logging while drilling tool, such antenna configurations can also be readily applied to wireline logging tools. Furthermore, the principle of reciprocity can be applied to obtain equivalent measurements while exchanging transmitter and receiver roles for each antenna. It is intended that the following claims be interpreted to embrace all such variations and modifications. 

1.-9. (canceled)
 10. An electromagnetic resistivity logging tool that comprises: a rotational position sensor; at least one receiver antenna; a first transmitter antenna having a first transmitter antenna orientation; and a second transmitter antenna having a second, different transmitter antenna orientation; and a processor to receive at least one of a phase and amplitude measurement in response to each of the transmitter antennas at each of a plurality of azimuthal angles to calculate at least one component of a coupling matrix using only transmitter antennas having the first and second transmitter antenna orientations.
 11. The tool of claim 10, wherein the processor calculates at least seven components of the coupling matrix using only transmitter antennas having the first and second transmitter antenna orientations.
 12. The tool of claim 11, wherein the processor calculates the seven components of the coupling matrix using only receiver antennas having the same orientation as said at least one receiver antenna.
 13. The tool of claim 10, further comprising a second receiver antenna, wherein said at least one of a phase and amplitude measurement is at least one of a phase difference between responses of the receiver antennas and an amplitude ratio between responses of the receiver antennas.
 14. The tool of claim 10, wherein the processor averages measurements in bins associated with each of the plurality of azimuthal angles and determines the coupling matrix components from the average measurements.
 15. The tool of claim 10, wherein the coupling matrix components are used to determine a log of a formation parameter for storage or display.
 16. The tool of claim 15, wherein the formation parameter is in a set consisting of anisotropy and vertical resistivity.
 17. The tool of claim 15, wherein the formation parameter is a formation dip angle.
 18. The tool of claim 10, wherein the processor determines a geosteering signal based at least in part on the C_(XZ) or C_(ZX) component of the coupling matrix.
 19. A logging method that comprises: receiving at least one of a set of amplitude and phase measurements as a function of azimuthal angle and tool position from a set of receivers having no more than two different receiver orientations in response to signals from a set of transmitters having no more than one transmitter orientation; determining from said set of measurements at least seven components of a coupling matrix as a function of tool position; finding at least one of vertical resistivity and anisotropy as a function of tool position from the coupling matrix components; and storing said vertical resistivity or anisotropy in the form of a log.
 20. The method of claim 19, wherein said at least one of a phase and amplitude measurement is at least one of a phase difference between responses of the receiver antennas and an amplitude ratio between responses of the receiver antennas.
 21. The method of claim 19, further comprising: averaging measurements in a plurality of bins associated with different azimuthal angles.
 22. The method of claim 19, wherein said determining is performed based on a linear system of equations that express a dependence of each receiver's response on azimuthal angle. 